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Abstract 

There exist magnetic fields in which particles bouncing between mirror 
points experience no net first-order guiding center drift. In such fields, 
even though the instantaneous gradient and curvature drifts are not zero, 
their total effect integrated over any bounce period vanishes, so that 

1 

particles merely wobble back and forth around fixed field lines. A class 
of two-dimensional drift-free fields, somewhat resembling the configu- 
ration found in the geomagnetic tail, is described; several proofs of the 
drift-free property are given, including some that suggests that the pro- 
perty of vanishing net drift might extend to non-adiabatic orbits. A general 
criterion for identifying drift-free fields is developed and a case of 
motion in a nearly drift-free field is also investigated. The theory is 
applied to the plasma sheet in the earth’s magnetotail and observational 
evidence is presented suggesting that the magnetic field there indeed 
approaches a drift- free configuration. 
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Introduction 


The purpose of this work is to examine a certain class of charged 
particle motions In a magnetic field* The motions involve "magnetic 
mirror" geometries in which two independent periodic inodes are possible 
gyration around field lines and "bounce" between mirror points - with 
two associated adiabatic invariants \x and J .The fields discussed 
here have the additional property that the net guiding center drift , 
averaged over any bounce period, vanishes; as will be seen, an analogous 
property can also exist for non- adiabatic orbits* 

An example of such a configuration is given in the next section: it 
has a finite curl and therefore it is not certain that particle distri- 
butions can be found which would provide the required current density , 
in the absence of any net charge transport. However, self-consistent 
plasma distribution functions certainly will exist for fields which 
are not drift -free but merely approach a drift -free configuration: in 
such cases the Integrated drift does not vanish but is merely small, so 
that the field can carry a substantial plasma density with only a modest 
current density# An analogous situation exists in dynamo theory where 
Cowling [i22L] proved the impossibility of axi symmetrical dynamos, 
y e ^ Braginskii £ 1964 a, b~J derived useful models by assuming a slight 
deviation from axial symmetry* 

A magnetic configuration which can be treated as a perturbed version 
of a drift-free geometry apparently exists in the plasma sheet of the 
earth's raagnetotail, explaining why the total current carried by this 
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sheet is relatively small, considering the amount of plasma in the sheet. 
The outer magnetosphere of Jupiter, as observed by Pioneer 10 \ Smith et 


al p , 1974 ] can also be represented locally by a two- dimensional drift -free 
field, to which a small component is added to preserve axial near- symmetry, 
so that it, too, qualifies as a nearly drift-free configuration. The field 


in Jupiter's outer magnetosphere might also be approximating an axi symm- 
etric drift -free geometry but, although such geometries are investigated 
in what follows, it is not known whether they actually exist. 


A Simple Example 


Consider a two-dimensional magnetic field of a form which is frequently 
used in models of the magnetotail plasma sheet 

B = B x (z) x + B z z (1) 

where 

B z - constant (2) 

Field lines of this configuration, in any plane of constant y and for 
the case where B x = >z, are given in Figure (l) , which also depicts 

(schematically) the projection of the orbit of a particle undergoing 
guiding center motion. One may now ask, does such an orbit undergo a slow 
drift in the y direction, orthogonally to the plane of the figure? 

This problem is easily handled by the mathematical device described 
below, although more direct proofs also exist 


and will be presented. 
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Suppose that an electric field 

5 “ e 0 2 ( 3 ) 

exists across the magnetic field:. (An electric field probably does exist 
across the plasma sheet, so that if equation (l) is assumed to represent 
a model of the magnetotail, equation ( 3 ) should probably also be assumed; 
this, however, is not relevant to the present calculation, in which the 
electric field is merely a mathematical device.) Let v be the bounce- 
averaged magnetic drift velocity of the given particle in the y direction., 
and let q be its electric charge; the particle then will increase its 
kinetic energy W at the rate q E Q v y per unit time. 

However, if the problem is transformed into a frame moving with the 
uniform velocity 

■a - £ VV M 

then in the non-relativistic approximation the electric field vanishes. 

If v y remains unchanged by the transformation (see below) then in the 
moving frame no energy gain is possible: this can only be reconciled 
with the previous assumption if v y « 0 , i. e . the particle’s guiding 
center "wobbles" around the guiding field line with no net drift. 

To see whether v y is conserved one should note that to the lowest 
relativistic order the magnetic field is not changed by the transformation 
to the moving frame and therefore any portions of v which depend 


on 
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B alone will also transform unchanged.* Ill addition v will contain some 

y 

terms involving the electric field ^ Northrop , 1963 , eq, 1,17 J . The 
above argument then states that v * 0 for any finite E and it is 

y 0 

plausible to assume (although this is not a rigorous proof) that the 
same also holds in the limit Eq ~ 0 , in which case the magnetic field 
is indeed drift-tfree* 

If this approach is accepted one obtains two additional insights 
into this type of motion* First, there is no need to assume that the motion 
is adiabatic . the same conclusions should also be valid for non** adiabatic 


motion, as long as v is suitably defined* An exception occurs in the 
limit B z ~ 0 , when (if) can no longer be used; indeed, in this limit mean- 
dering modes of motion exist ^ Sonne r up , lgjlj which have finite v y . 


Nearly Drift-Free Fields 


A second advantage of the above approach is that it allows one to estimate 
V y some cases i* 1 which the geometry is only approximately drift -free* 

Suppose B a is not constant but varies slowly with x ; then the 
preceding argument no longer holds, since u is no longer a uniform 
velocity. Assume now 

V 0) = 0 (5) 

(as in the field in Fig. l) and consider the adiabatic motion of a charged 
particle in the plane z - 0 « The electric drift will carry such a 
particle in the x direction and due to conservation of the magnetic 
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moment p the particle f s kinetic energy W will change -with x at 
the rate 

= n x (6) 

z 


By previous arguments the rate with which W changes in time equals 

q E Q v^ and it can also he found by multiplying (6) with the electric 
drift velocity , giving 

v y - (n/qB^) TIB z /Qx (7) 

This can also be written as 


v 

y 




( 8 ) 


vhere v is the particle’s velocity, R its gyration radius and L 

6 x 

the scale length of the variation of B z . By way of comparison, the 

instantaneous guiding center drift velocity is of the order of v R /l 

g * 

where L is the much smaller scale of variation of the total field. 


The derivation of v^ for adiabatic particles not confined to 
the equatorial plane is much more difficult and requires the bounce- 
averaged Hamiltonian p, j), discussed later on. 


Alternative Proofs 


Consider Newton's equations of motion for a particle of charge q 
and mass m in the field of (l) 0 The x component of the acceleration 


is 
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x = 1 B z y / m (9) 

and It integrates to 

x = q B z (y - y Q )/m (lo) 

If the motion in x is periodic, as happens in the adiabatic case, 
y will also vary in a periodic fashion. Even if the motion is non- 
adiabatic but x varies in a periodic or almost periodic fashion, 
the same property will hold true for y and no net drift will occur, 
showing that the motion is drift-free. As before, the argument cannot 
be used in the limit - 0 , since x vanishes in that case, and as 
has already been noted, the motion in that limit is indeed not drift- 
free. 


A more conventional proof of the drift-free character of adiabatic 
motion in the field (l) involves the averaging of the instantaneous 
guiding center drift over a bounce period. The calculation is lengthy 
and is therefore given in the appendix; it is assumed there that the 
magnetic configuration possesses mirror symmetry with respect to the 
plane z = 0 but, as is shown in what follows, this requirement is 
not essential. 

The Averaged Hamiltonian 

A powerful tool for treating motions in a magnetic mirror geometry is 
provided by the bounce-averaged Hamiltonian K( o<, ^ j), introduced 

by Northrop and Teller [aaffi] . Its use represents a second averaging , 
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over the periodicity of the mirroring motion, in addition to the 
gyration averaging which leads to the guiding center theory* 


Let he Euler potentials of the magnetic field (e.g. 

review by Stern [mo] ) and let a particle move in this field and 
conserve the two lowest adiabatic invariants \x and J * In the absence 
. of electric fields the Hamiltonian K(o(,p, p,, j) for the averaged 
motion then equals the particle f s kinetic energy W „ 


Attention will be confined to cases in which, due to some symmetry, 

K does not depend on - in particular, either to 2- dimensional 
fields where 

B = B(x, z) ; B^= 0 ; f - y 

or to axially symmetric fields where (using cylindrical coordinates $ , 
z and y ) 

B = B( f, z) ; = 0 ; p = 'f 

The drift velocity (combining gradient and curvature drifts) then 
parallels the direction of and its average value can be found by one 

of Hamilton's equations of motion 

» 

< P > = (c/e) k/ (11) 


In particular, if K is independent of both and p , the average 
drift in the direction of vanishes and the geometry will be drift- 
free. This is a general condition for drift-free geometries but it does 
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constitute a useful criterion for identifying such fields, since the 
explicit form of K is difficult to derive [e„g. Chen and Stern . 1975 j 
We therefore proceed to develop an alternative criterion. 


The Euler potentials are related to B through 

B = Vo(* Vf ( 12 ) 

Let it be first assumed (an assumption which will later be removed) 
that c( has mirror symmetry with respect to the plane z « 0 , i.e. 
for a 2- dimensional field 

z) = UU,-z) (13-a) 

and for an axisymmetric field 

oU 2)= -2) (13-b ) 


If a particle guided by the field line ( o( 

kinetic energy tt^Yadiabatic invariants (pi, j) and^ reflected at 

z = z , then 
m 

J = (4/m) (2W) 1 / 2 

jD y. / y 1 

m last equation implies that a function f(z o() 

nr 

exists such that 


I 


m 


(1 - 


B/fe ) 
m 


1/2 


( Be/ dz)| dz (l4) 


, p ) of this field has 


J W 


-1/2 


= f (z , oC ) 


nr 


( 15 ) 


The magnetic moment 


M- t on the other hand, satisfies 
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w/n = B m = B (z m , o< ) (16) 

Elimination of z^ between the last two equations implicitely gives 
the functional form of W( 4, n, j) which, in the absence of electric 
fields, equals the Hamiltonian K of (ll)» If W is independent of 

q 4 , this means that such an elimination removes oi as well, which 
will only happen if the right hand expressions depend on each other, 
i.e„ 

f(z m , <*,) - f(b( V o<)) (IT) 

This is easier to express if (11*) is transformed to a form in which 
the field intensity B becomes the integration variable, namely 
B 

l(c ^ V = y (1 “ B / B m }l/2 ^s/^DBj dB = G(B m ) (18) 

Here = B(z=0,cO is the equatorial field intensity on the field 

line with the given value of oC and the last equality, involving some 

function G independent of o( } must hold if the field is to be drift - 
free* 

A simple extension of the arguments of the section dealing with 
nearly drift-free fields shows that B eq may not depend on U because 
if such a dependence existed, then in the limit J — ► 0 a particle 
would experience a net drift, In that limit a particle is confined to 
the plane z = 0 : assuming an electric field 

" E o^ 


E 


(19) 



-li- 
lt is found that the electric drift of such a particle would he in the 

direction of • If B were allowed to vary with then the 

eq 

particle's energy would change in the course of its electric drift due 
to the conversation of \i and, as outlined earlier, this would imply 
a net magnetic drift in the direction of « 

The condition that I is independent of o< thus reduces to the 
requirement that IIs/DB at constant cC depends on B alone - i.e. 

s « s(b) (20) 

Another way of expressing this relationship would he in the form of 
a Jacobian 

Tl(s, B)/Q(q i , q^) - 0 (2l) 

where (q , q^) equals on (x, z) or ( S , z) depending on the geometry. 

In the cases discussed here both s and B can he expressed in terms 
of the Euler potential o( and its derivatives and fcl) then reduces 
to a rather complicated integro-differential equation. 

In all this the requirement of symmetry expressed by equations (13) 
v , provided (3 ) holds j 

can be waivedJT The reasoning behind this is outlined below: it can he 
readily expressed in conventional mathematical terms hut doing so has 
little real effect beyond lengthening the presentation. 

Suppose that the symmetry condition holds and consider two points 
P 2 = (x, z) and « (x, -z) in a two-dimensional field, symmetrically 
located on the same field line at equal distances from the plane z - 0 
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(for axial symmetry replace x everywhere by $ ). One now readily 
finds (e.g* from equations (a- 1) and (a- 2) in the appendix) that both 
the curvature drift and the gradient drift at P and at P g are the 
same. By extension it can also be shown that the average drift of either 
kind experienced by the particle along the portion of its guiding field 
line above z = 0 is the same as the corresponding average drift below 
Z = 0 * S±uce the ^erage gradient drift in any period cancels the average 
curvature drift, it also follows that a similar canceUat ion holds just 
for the portion of the orbit above z « 0 and just for the portion below 

2 ~ 0 ■ Thus * lf one takes different field line shapes satisfying (l), 
(2) and (5), and splices them together at z = 0 , the motion in the 
resulting field will still be drift-free, as has been claimed. 

In the case of the field of (l) and (2) it is possible to choose 

z 

' = x B z - J B x dz (22) 

0 

f = y (23) 

calculation then shown that both B and s depend on z alone, the 
elimination of which leads to the condition (20). Analytic examples or 
even existence proofs for the axially symmetric case are not known. 

Application to the Plas m a Sheet 

The form of the magnetic field given in equation (l) and in particular 
the configuration shown in Figure 1 resemble in many way the magnetic 
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field observed in the earth's tail, in the plasma sheet region* In 

that region the value of B gradually decreases from about 4-7 at 
^ z 

20 R e to 0*5 7 at 70 Rg ^ Behannon , 1970 * table on p* 751 ; Boiling 
and Wolf , 1214] , while « 4 7 /R e at a distance of 30 R g • 

The preceding theory suggests that the mean velocity v of charged 

y 

particles across the plasma sheet is relatively small* If the current 
density observed in the tall can be entirely attributed to the mean drift 
of particles « a point examined separately later on - then the magnitude 
of the tail current can be used to estimate Vy in the following way. 

At a distance of 20 R the change in B across the plasma sheet 

5 

amounts to about 20 7 , representing a current of about 10 amperes for 

each R 0 of plasma sheet length. Assuming that the current is carried 
mainly by protons (which are considerably more energetic than electrons 
and therefore drift more rapidly) and taking the plasma density in the 
middle of the sheet as 0.5 protons /cm 3 over a width of 6 R e , one 
finds that in order to maintain the observed current, the protons must 
move with an average cross-sheet velocity 

v « 5*13 km/sec 

y 

A very rough theoretical estimate of what v ought to be is 

y 

obtained from eq* (8) , which assumes protons to be adiabatic and 
equatorial* If a typical plasma sheet proton has an energy of 2 kev, 
one finds v « 620 km/sec, R g « 1/3 R e (in a field of 57 ) while the 
work of Behannon [ms] suggests 1^ « 15 R e 3 leading to 


Vy » 6*9 km/sec 
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The agreement obtained here should only be considered as qualitative, 
for at least two reasons. First, the conditions for which (8) has been 
derived are not generally met and secondly, we have ignored the effects 
dtte to density variations in a diamagnetic plasma. 

Concerning the second point it should be noted that the plasma sheet 
current density is contributed not only by the drift velocity v 

y 

of charged particles but also by the variation in the density M of 
magnetic moment. The latter contribution is called the magnetization 
current density £e,g. Longmire , iss] and equals V * M • 

In a strictly drift-free configuration of the type discussed here, 
if the plasma density N does not vary with x then V* M does 
not contribute to the total current, for the following reason . Assume 
that the plas m a sheet is bounded at z = t Zq by plasma- free regions 
(Figure 2). If V XM is now integrated over a rectangular cross- 
sectional area of the sheet, bounded by z « t (z > z^), x — x 
and x « x g , then the integral equals the total magnetization current 
through the given area. Converting the integral into a contour integral 
of M around the area, it is seen that as long as no x dependence 
exists, the result vanishes. 

Even if N does not depend on x , however, a weak dependence of 
on x (as was assumed here) causes a non-vanishing total magneti- 
zation current, since the contributions of the segments with x = x 
x — X 2 ^he contour integral then no longer cancel each other® 

To the lowest order of approximation the result turns out to be propor- 
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tional to ; in a similar way, the current conducted "by 

particle drifts is also proportional to /^»x , which suggests 

that in nearly drift- free geometries the contributions of both these 
mechanisms to J are roughly of the same order. 

It is interesting to compare these results to the work of Bird 
andBeaxd [l^] who computed numerically the contributions to j 
from the curvature drift current, the gradient drift current and the 
magnetization current, in a model field representing the magnetospheric 
tail (Figure 1 in their paper). As expected the first two sources of 
current tend to have opposite effects; the magnetization current has 
an appreciable effect on the dependence of J. on z , but its 

contribution to the total current is small and the total current density 
from all sources is likewise small. 

One also should be careful here to note (as was pointed out by one 
of the referees) that the guiding center approximation cannot always 
be assumed to hold. The mean energy of plasma sheet protons near z « 0 
at 18 R e peaks between 2 and 5 kev [game et a 1., 1967 ] . Taking 
5 kev as a typical value of the energy at 30 R e where B z « 27 , we 
find a gyration radius of about 0.6 R 0 . On the other hand, if 
^ B x^ 2 ~ scale of field variation is about 

B z / (^B x /^z) » 0.5 R e 

This suggests that near the mid-plane of the plasma sheet many of the 
protons may be moving in non- adiabatic orbits. The preceding theory 
suggests that the motion of such protons, too, is nearly drift-free, but 
for more quantitative conclusions their orbits must be examined numeri- 
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caTly, as has been done in the work of Pudovkin and 
Numerical work can also clarify the transition from the drift -free mode 

to the meandering mode of Sonnerup [l97l] which, as was noted, is not 
drift free • 
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Appendix 

For completeness the vanishing of the mean drift velocity in the field 
of eq. (l) is derived here by straightforward integration. The gradient 
and curvature drifts are 


V SS 

-g 

(v A 2 /2WB 2 ) B x 7B « - (v i 2 /2<y)(B 2 /B 3 ) dB /dz y 

(A-l) 

V = 

—c 

- (v^ 2 /ttB 2 ) B X [b x ( 7 x B)J 


= 

(v a 3 /6)) (B 2 /B 3 ) dB^dz £ 

(A-2) 

Assuming that the particle oscillates between mirror points 

with a 

period 

T , the mean drift velocity is 



<y> = (4/r) j (y) dt 

* av 

(A-3) 

where the integration extends from z « 0 to the mirror point and ( ) 

av 

implies averaging over a gyration period. It will be found convenient 

to use 

as integration variable the field intensity B observed by the 

particle and to express eq, (A-3) as 



<y> = (Vf) f m (y)„_ (dB/dt)- 1 dB 

B av 

z 

(A-k) 

To evaluate dB/dt two unit vectors are introduced, B aligned with 

B and 

— orthogonal to in in the (x, z) pi»r, 0 . 



B = B -1 (B x + B z) 

X z — 

(a-5) 


- = B_1 ^" B z X +B X z) 

(A-6) 
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This gives 


A 

z 


B* 1 (B B + B N) 

z — X — 


(A- 7 ) 


Using the z coordinate to measure the particle r s progress along its 
guiding field line one gets 


dB/dt 
“by (A- 7 ) 


(dB/dz)(dz/dt) = (B /b) dB /dz (v*z) 

X X — 


(a- 8) 


■A 

v*z 


B* 1 (B v + B v ) 
z ^ x n 


(A-9) 


vhere v q is the velocity component along M . Averaged over a guiding 
center gyration, to lowest order 


Cv ) =0 

n av 


(A-10) 


and therefore 


(dB/dt) 


av 


Cv '/ B * B z / b2) 


(A-ll) 


Also, if the rapid gyration around the guiding center is averaged out, 
then (y ) w of equation (a- 4) simply becomes the combined drift velocity 
contributed by (A-l) and (a-2) 


(y) av - dB x /dz (20>B 3 r 1 (2 B* vf - B| ) 


(A- 12 ) 


Denoting by the field intensity at the mirror point, one finds 

from the conservation of the magnetic moment 
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vf - v 2 (B/B m ) 

f = v“ [!■- (B/B m )J 


(a- 13 ) 


Finally, if (A-U) to (A-13) are substituted in (A-4), the resulting 
expression has the form 


<;,y> = const, 


f B m 2 B 2 B - B 2 B - B 3 
/ z m z 

Bz (b 2 - b 2 ) 1 ^ 2 (b - b) 1 / 2 


M 

B 2 


(A-3A) 


The constant appearing here contains T, u)/b , v and B 5 the 
~ z 

derivative dB /dz cancels out while B itself is expressed in 

X X 

terms of B , as in the first factor in the denominator. It may then 
he shown that the integral vanishes, hy dividing it into two parts and 
using integration by parts to show that they cancel. It follows that 
the mean guiding center drift vanishes, as has been claimed. 
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Captions to Figures 


Figure 1 - Field lines in a drift-free magnetic field, with the 

projection of a trapped orbit schematically- shown 0 

Figure 2 - Integration contour used in showing that in the absence 

of an x-dependence of plasma density and magnetic field, 
the contribution of the magnetization current to the 
total current in a drift-free model of the plasma sheet 


vanishes 
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